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I. INTRODUCTION 

An important question in present-day cosmology con- 
cerns the expansion history of our universe. Recent mea- 
surements indicating a late-time period of accelerated ex- 
pansion [f] have been met with great interest. Possible 
theoretical explanations for this observation typically in- 
volve new fundamental physics including a cosmological 
constant [2] , quintessence-type models [3] with one [4] or 
two [5] scalar fields, k-essence [6], and exotic equations 
of state like that of the generalized Chaplygin gas [7] . 

Other astrophysical observations claim evidence for a 
time-dependent fine-structure parameter [8] . Early spec- 
ulations in the subject of varying couplings date back 
to Dirac's large-number hypothesis [9]. Subsequent the- 
oretical investigations have shown that a spacetime de- 
pendence of the fine-structure parameter arises naturally 
in candidate fundamental theories and is often accompa- 
nied by variations of other gauge or Yukawa couplings 
[10]. In light of these facts, a confirmation of the exper- 
imental observations and the search for realistic models 
that permit other particle-physics and cosmological pre- 
dictions have assumed particular urgency [If]. In this 
context, studies along the lines of the Bekenstein model 
[12] and its generalizations have suggested dark matter 
and a cosmological constant [13], an ultra-light scalar 
field [f4], and quintessence [15] as driving entities for a 
varying fine-structure parameter. 

The presence of varying couplings implies a breaking of 
invariance under temporal or spatial translations. This 
can be seen as a special case of the violation of space- 
time symmetries, which also include Lorentz and CPT 
invariance. In fact, time-dependent couplings typically 
affect these additional symmtries as well [16]. Note also 
that Lorentz and CPT breakdown has been suggested in 
a variety of approaches to fundamental physics. We men- 
tion string theory [17], spacetime foam [18,19], nontrivial 
spacetime topology [20], loop quantum gravity [21], and 
noncommutative geometry [22]. Lorentz and CPT vio- 
lation provides therefore another independent signature 
for an underlying theory. 

Scalar fields are a common feature in all of the above 
contexts. This shows that scalars play a key role in the 
search for fundamental physics. It becomes therefore 
interesting to investigate cosmological expansion, vary- 



ing couplings, and Lorentz violation within a single can- 
didate fundamental model containing scalars. In the 
present work, we shall consider an iV = 4 supergravity 
model in four spacetime dimension. The model contains 
two scalar fields, an axion and a dilaton, with a potential 
that we model with mass-type terms. This framework, al- 
though not fully realistic in its details, incorporates many 
features expected to be present in an encompassing the- 
ory: the pure four-dimensional iV = 4 supergravity is a 
limit of the iV = 1 supergravity in 11 dimensions, which 
is contained in M theory. 

The paper is organized as follows. Section II describes 
the basics of our model. In Sec. Ill, we demonstrate 
how, for a suitable range of parameters, a late-time pe- 
riod of accelerated cosmological expansion can arise in 
such a supergravity model. The associated variation of 
the fine-structure parameter and the electromagnetic 9 
angle is investigated in Sec. IV. Section V discusses the 
violation of Lorentz and CPT symmetry in the presence 
of spacetime-dependent scalars. A brief summary is con- 
tained in Sec. VI. 



II. BASICS 

The starting point for our investigation is an isotropic 
homogeneous flat {k = 0) Friedmann-Robertson- Walker 
universe with the usual line element 

ds^ ^dt^ ~a^{t) {dx^ + dy^ + dz^) , (1) 

where a{t) denotes the scale factor and t the comoving 
time. The assumption of flatness simplifies the analysis 
and seems justified in light of recent measurements [23]. 

In the present context, we are interested in physical 
effects past the formation of Hydrogen. We therefore ne- 
glect the energy-momentum contribution of radiation to 
the source of the Einstein equations. Galaxies and other 
matter are modeled in a standard way by the energy- 
momentum tensor of dust: 

Tf,i, = pui_,u^ . (2) 

Here, p is the energy density of the matter and is a 
unit timelike vector orthogonal to the spatial hypersur- 
faces, as usual. 
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In addition, our model contains two real scalar fields, 
an axion A and a dilaton B, coupled to the electro- 
magnetic field F^^. The fields are described by the La- 
grangian density £sg obeying 



■sg 



where 



-\k^MF^,F^''' - Ik^NF^^F^''' , (3) 



M = 



N = 



(1 + + £2)2 _ 4^2 

A(A^ + B^-1) 

(1 + A2 + B2-)2 _ 4^2 



(4) 



Here, F^^ = ei"'P''Fp„/2 denotes the dual field-strength 
tensor and g = — det(5^,^). For convenience, we rescale 
p p/k and F^"" F^""/^ in what follows. Then, 
the gravitational coupling k docs not appear explicitly in 
the equations of motion. Under the assumption that the 
dust arises through fermions uncoupled from the scalars 
A and B and with the identification of F^^" with one of 
the graviphotons, this model fits into the framework of 
N = A supergravity in four dimensions [24,16]. 

The pure N=4 supergravity is known to be incomplete. 
For example, its spectrum differs from the observed parti- 
cle species in nature, and the matter couplings of F^" are 
nonminimal. This latter issue can be avoided by gauging 
the internal SO (4) symmetry leading to a potential for 
the scalars that is unbounded from below [25]. However, 
besides additional fields, any realistic situation requires 
a stable vacuum. We are therefore led to a phenomeno- 
logical approach and take the effective potential to be 
bounded from below. We further assume that the scalars 
remain effectively uncoupled from the matter. This may 
explicitly break the supersymmetry, but does not repre- 
sent a problem phenomenologically. 

To lowest order, we model the potentials for these 
scalars by quadratic self-interactions for the axion and 
dilaton parametrized by the coefficients niA and ms, re- 
spectively. Then, our full Lagrangian density £ reads: 

£ = - l^{m\A^ + mlB^) + g^^T'"' . (5) 

Note that ruA and ms cannot be identified directly with 
the masses of A and B because of the non-canonical ki- 
netic terms. 

The next step is to determine the equations governing 
the evolution of the universe. Our above assumptions of 
matter dominance, homogeneity, and isotropy imply that 
we can neglect electromagnetic effects for the large-scale 
dynamics of the universe. We therefore set F^,^ — for 
the moment. Moreover, the matter density p and the 
scalars A and B are (like the scale factor a) functions of 
the comoving time t only. Then, the respective 00 and jj 
components of the Einstein equations including the dust 
are given by: 



m\A'^ + tubB'^ + 



A-+2^ = m\A^ + mlB^- 



A^ + B^ 

2E2 
A^ + B^ 

2S2 



+ 2p 



(6) 



where the dot denotes a derivative with respect to the 
comoving time. The off-diagonal components arc trivial. 

The equations for the time evolution of the axion and 
the dilaton are obtained by variation of C with respect 
to A and B, as usual: 



= -2a^m\A 




-2a'mlB-^iA'+B'). (7) 



The above equations of motion for A and B imply that 
the energy contained in the scalars is conserved sepa- 
rately: 



= a^D^e^o 



d_ 
di 



a 

4B2 



(i2 + 

1 



-{m\A^ + mlB^) 



2^2 



(i2 ^ ^2^ 



m\A 



2 m%B^ 



(8) 



where O^'^ denotes the symmetric energy-momentum 
tensor associated with the scalars A and B. This es- 
sentially arises as a result of the fact that our model fails 
to contain interactions of the axion and the dilaton with 
the fermion dust. This is a phenomenologically reason- 
able assumption: energy exchange between the scalars 
and the dust can only be mediated by the electromag- 
netic field. However, the matter content of the universe 
is electrically neutral and stable against photodecay on 
macroscopic scales, which excludes the required rates of 
photon exchange. 

As an immediate consequence of the above result, the 
energy-momentum tensor T^^ of the dust is also covari- 
antly conserved. This gives the equation d{pa^)/dt = 0, 
which can be integrated to yield p{t) = Cn/a^(t). Here, 
the integration constant Cn = PnO,^ is determined by the 
matter density pn = p{tn) and scale size Un = a{tn) of 
the universe at the present time i^. 

Note that Eqs. (6) and (7) are four equations for the 
three miknown functions a{t), A{t), and B{t). To demon- 
strate the dependency among them, one can proceed as 
follows. Multiplication of the first Einstein equation in 
(6) with and subsequent application of a time deriva- 
tive eliminates the term containing p by energy conserva- 
tion. Moreover, with the aid of the equations of motion 
for the scalars in the form of Eq. (8), one can rearrange 
the remaining expression to obtain the second one of the 
Einstein equations in (6). 
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III. ACCELERATED COSMOLOGICAL 
EXPANSION 

In this section, we investigate the sohitions of the Eqs. 
(6) and (7), which govern the dynamics of our model. 
The two Einstein equations (6) imply 



6^ 

a 



~B2 



- P ■ 



(9) 



Thus, in a physical model at least one one of the param- 
eters niA or TUB must be nonzero to yield an accelerated 
expansion rate a{t) > 0, as expected. This is also con- 
sistent with a previous analysis taking uia = tub = 0: 
in this case, all equations of motion can be integrated 
analytically leading to the conventional decelerated ex- 
pansion ^ t^/^ of a matter-dominated universe in the 
present cosmological epoch [16]. 

In special situations in which niA, or ms, or both are 
nonzero, some analytical results can be obtained. These 
results together with some general comments on the spec- 
trum of solutions can be found in Appendix A. However, 
for the present purposes numerical methods appear most 
promising for further progress. A phcnomcnologically in- 
teresting input-parameter set should minimize the devi- 
ation from experimental data. As a result of the nonlin- 
earity and the complexity of the present system of dif- 
ferential equations, it appears difhcult to find a global 
minimum. A numerical search has given a variety of pa- 
rameter sets optimizing locally the departure from obser- 
vations. One such set sufficient for our purposes is the 
following: 

2.7688 X IQ-''^ GeV , 
3.9765 X IQ-'^^ GeV , 
2.2790 X lO-^"* GeV^ , 
1 , 

1.0220426 , 

-8.06401 X 10"^^ GeV , 
0.016598 , 

-2.89477 X 10"''''^ GeV . 



rriA 

THb 

Cn 

a{tn) 

i(tn) 
B(tn) 
B(in) 



(10) 



The numerical analysis also showed that there are two 
sets of solutions: one in which the scale factor a{t) in- 
creases with time and another one in which the scale 
factor decreases. In what follows, we only consider the 
observationally relevant solution with a{t) > 0. 

Our next step is to compare the solution determined 
by the parameters (10) with experimental data. Mea- 
surements show that the expansion of our universe is 
consistent with a flat "canonical" model containing non- 
relativistic matter (dust) and a cosmological constant A. 
The parameters for this canonical model arc inferred from 
observations of high-redshift supernovae [26] : 

Qm = 0.30 ± 0.04 , 
Oa = 0.70 ± 0.04 , 

ilo = (70±4)km-s-^-Mpc-^ . (11) 



Here, and are the respective energy densities as- 
sociated with the matter and the A term relative to the 
critical density 3Hq^/k. The Hubble constant is denoted 
by Hq, as usual. To make contact with observations, we 
can therefore compare our solution determined by the in- 
put values (10) with the above canonical model (11). To 
allow only small deviations of the supergravity cosmol- 
ogy from the best-fit canonical model, we have chosen the 
variations in the parameters (11) to be somewhat smaller 
than the experimental errors [26] . We also take the vari- 
ations of the relative energy densities Om and Q\ to be 
constrained by f^M +f^A = 1, so that the canonical model 
yields a fiat (k = 0) universe for all parameters (11). 

The supergravity cosmology leads to a value for the 
Hubble constant of Hq — 69km-s~^-Mpc~^ consistent 
with observations. The canonical model yields tn = 
1.35 X 10^" yr for the present age of the universe. Our 
supergravity cosmology implies a similar age of tn = 
1.30 X lO^Oyr. 
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FIG. 1. Time evolution of the scale factor a{t) and its sec- 
ond derivative d{t). The solid and dashed lines correspond 
to our supergravity universe and the canonical model, respec- 
tively. Note that for approximately the second half of its 
lifetime, the expansion of the universe is speeding up in both 
models. 



Other features of the cosmological expansion are illus- 
trated in Fig. 1. The upper panel shows the evolution 
of scale factor from the Big Bang to the present. The 
time dependence of the scale-factor acceleration a{t) is 
depicted in the lower panel. The solid and the dashed 
lines refer to our supergravity cosmology and the canon- 
ical model, respectively. The fractional comoving time is 
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determined by the corresponding values of each model 
given in the previous paragraph. Note the recent accel- 
erated expansion of the universe in both cosmologies at 
comparable rates. 

Although the differences in scale-factor evolution bet- 
ween our and the canonical model seem relatively mi- 
nor in Fig. 1, it is necessary to investigate whether 
these deviations are such that consistency with obser- 
vations is maintained. Experimentally, the time depen- 
dence of the scale factor is inferred from luminosity and 
redshift measurements of distant cosmological objects. 
Conventionally, the luminosity observations are trans- 
lated into a luminosity distance that can be ex- 
pressed as a distance- modulus value (m — M). In Fig. 
2, A(m - M)n=o = (m - M) - (m - M)n=o is plot- 
ted as a function of redshift z. To emphasize the accel- 
eration effects, we have subtracted the distance modu- 
lus (m — M)q^q of an empty (k — —1) non-accelerated 
universe with a{t) = Hot and Ho = 70 km-s~^ -Mpc"^. 
As before, the solid line corresponds to our supergravity 
model and the dashed line to the canonical model. The 
dotted line represents the empty universe. The shaded 
area corresponds to the parameter space (11). In the 
shown redshift region, which corresponds to the current 
observational range of z < 1, both models explain the 
supernovae data. At higher rcdshifts z > 1, the models 
could in principle be distinguished by future experiments 
of this type. 
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FIG. 2. Distance modulus relative to an empty universe 
A(m — M)n=o versus redshift z. Our supergravity cosmology 
is represented by the solid line and the canonical model by 
the dashed line. The dotted line corresponds to the empty 
universe. The shaded region marks the parameter range (11). 

Next, we briefly discuss the time evolution of the 
scalars A{t) and B{t). The values (10) lead to the 
comoving-time dependence of the axion and dilaton fields 
depicted in Fig. 3. For clarity, the value of B{t) has been 
multiplied by a factor of ten. The scalars A{t) and B{t) 
vary significantly only in the early universe. At later 
times, they display a comparatively constant behavior, 
which is phenomenologically necessary because the ax- 
ion and the dilaton determine the variation of the fine- 
structure parameter, which is tightly constrained experi- 
mentally. A more detailed discussion of this topic can be 



found in the next section. 

We conclude this section with a few comments on 
two further quantities of interest. One of them is the 
equation-of-state parameter w = Ps/Psj where denotes 
the pressure of the scalars and Ps their energy density. 
Figure 4 shows a plot of w versus the fractional comov- 
ing time. Note that w is significantly different from — 1 
only at early times. It follows that in the recent past the 
scalars affect the expansion of the universe essentially in 
the same way as a cosmological constant. 
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FIG. 3. Time evolution of the scalars A{t) and B{t) at 
early cosmological times. In the recent past of our model uni- 
verse, which is not shown here, A{t) and B{t) are essentially 
constant. 
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FIG. 4. Time evolution of the equation-of state parameter 
w. At late times, w ^ —1, so that the scalars essentially obey 
the cosmological-constant equation of state. 

The second quantity of interest is the distribution 
of energy between the nonrelativistic matter and the 
scalars. In Fig. 5, the relative energy density of the dust, 
r^M, is plotted versus the fractional comoving time. The 
relative energy density Qs associated with the scalars 
obeys J^m -I- = 1 in a flat universe, so that the shaded 
region corresponds to the time evolution of ilg. At late 
times, ris becomes the dominant energy component in 
our model. This behavior, which is characteristic for a 
conventional universe with positive A energy, arises be- 
cause the equation of state for the scalars approaches 
that of a cosmological constant. The supergravity cos- 
mology exhibits an energy distribution differing from the 
conventional case only at early times: f2s also dominates 
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close to the Big Bang. However, we remind the reader 
that our simple model neglects important effects prior to 
the recombination time. 
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FIG. 5. Relative energy density of matter Qm versus frac- 
tional comoving time. The shaded area shows Qs, which 
corresponds to the energy associated with the axion-dilaton 
background. At late times, Qs dominates, which parallels the 
cosmological-constant situation. 



IV. VARYING COUPLINGS 

In this section, we consider excitations of F^^, in the 
axion-dilaton cosmology discussed above. In most ex- 
perimental investigations, the spacetime regions involved 
are small on a cosmological scale, so that it is appropri- 
ate to work in local inertial frames. A few of the fol- 
lowing results have been derived previously for the case 
TTiA = niB = [16]. However, we discuss them here 
anew, both for completeness and to emphasize qualita- 
tive differences to the present case. 

The conventional electrodynamics Lagrangian density 
in a local inertial frame can be taken as 

1 f) 
/[e'^^'' 167r2 '"^ ■ ^ ' 

We remind the reader that the 6'-angle term fails to con- 
tribute to the classical Maxwell equations as long as 9 
remains constant. Comparison with our supergravity 
model yields = 1/M and 6 = Att'^N. Note that M 
and N are functions of the scalars A and B. It follows 
that the time dependence of A and B in the supergravity 
cosmology induces a time variation of e and 9. Thus, in 
the present model the fine-structure parameter and the 
electromagnetic 9 angle acquire related spacetime depen- 
dences in a general coordinate system. 

In Fig. 6, the relative variation {a — an)/an of the 
fine-structure parameter a = e^/47r is plotted versus the 
fractional look-back time 1 — t/tn to the Big Bang. Here, 
a,! — 1/137 denotes the present value of the electromag- 
netic coupling. The solid line represents our supergravity 
model with the input (10). Also plotted is the recently 
reported Webb dataset [8] obtained from observations of 
high-redshift spectra. Both the supergravity cosmology 
and the data show a relative variation of a by parts in 



10^. Although our model fails to provide a good match 
to the data, it does exhibit nonlinearities in a{z), a nec- 
essary feature to fit all experimental constraints. 
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FIG. 6. Relative time variation of the electromagnetic cou- 
pling versus fractional look-back time to the Big Bang. The 
present model with the parameters (10) is represented by the 
solid line. Also shown is the Webb dataset. 

We mention that a generic choice of model parameters 
typically yields a much larger variation of a, which can 
be understood as follows. The supergravity cosmology 
is governed by the coupled equations of motion for the 
scale factor a{t) and the scalars A{t) and B{t). The ob- 
served value of the Hubble constant provides the exper- 
imental constraint d/a 10"^*^ yr~^. One would there- 
fore expect that this is also the approximate variation 
of the axion-dilaton background. This value, however, is 
roughly six orders of magnitude larger than the scale for 
a/ a suggested by the Webb data [8]. 

In a general phenomenological model involving some 
scalar field (j), the difference in the scales for d/a and 
a/ a is typically bridged as follows. One normally ex- 
pands a(0) = an{l + $0) + . . ., which permits adjusting 
the parameter ^ to obtain the desired order of magnitude 
for a /a. Note that this normally introduces an additional 
scale into such a phenomenological model. In the present 
supergravity framework, however, the dependence of a 
on the scalars A and B is fixed, so that adjustable coef- 
ficients in the electromagnetic coupling are absent. The 
above discussion suggests that this is a critical issue in 
models with parameter-free interactions linking cosmo- 
logical expansion and variation of couplings via the same 
scalar field(s). 

V. APPARENT LORENTZ VIOLATION 

It is known that spacetime-varying couplings are typi- 
cally associated with Lorentz violation [16]. We summa- 
rize this result first before further investigation. In the 
present context, Lorentz violation is exemplified by the 
modified Maxwell equations resulting from the effective 
Lagrangian density (12): 

- + A(a^^)i^p. = . (13) 
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For spacctimc-indcpcndcnt e and 9, the usual electro- 
dynamics equations are recovered. In our axion-dilaton 
cosmology, however, the gradients of e and appearing 
Eq. (13) are nonzero, approximately constant in local 
inertial coordinates, and act as a nondynamical exter- 
nal background. This vectorial background determines 
a preferred direction in the local inertial frame violating 
particle Lorentz symmetry, as defined in Refs. [30,39]. 

We remark that general Lorentz and CPT violations of 
this type are described by the Standard- Model Extension 
(SME) [30-32]. This effective-field-theory framework 
contains all coordinate-independent Lagrangian terms 
formed by contracting Standard-Model operators with 
Lorentz-violating parameters. For example, the term in 
Eq. (13) involving can be identified with the fc^^ op- 
erator in the SME, which breaks both Lorentz and CPT 
invariance. We also mention that the SME provides al- 
ternative explanations for the baryon asymmetry in our 
universe [33] and for the observed neutrino oscillations 
[34]. 

Particle Lorentz violation through couplings varying 
on cosmological scales is independent of the chosen refer- 
ence frame and is not only a feature in electrodynamics. 
Consider general equations of motions in a given space- 
time region. If they contain the nonzero gradient of a cou- 
pling in a particular set of local inertial coordinates, the 
gradient is present in any coordinate system associated 
with the spacetime region in rjucstion. Although dynam- 
ical quantum scalar fields are theoretically attractive as 
driving entities for the variation of couplings, the above 
argument for Lorentz violation is also valid for classical 
scalar fields. In fact, the coupling need not be associated 
with a dynamical field at all. 

In what follows, we shift our focus from the electromag- 
netic sector to the time-dependent scalars themselves. In 
an effective vacuum characterized by sizable cosmologi- 
cal solutions A{t) and B{t), the propagation of the ax- 
ion and the dilaton is modified relative to the constant- 
background situation. To investigate local effects of the 
background solutions on the axion and the dilaton at a 
certain spacetime point xq one can proceed in a local 
inertial frame. Propagating disturbances in the axion- 
dilaton background are described by small perturbations 
SA and 6B away from the cosmological solution: 

Aix) Aix) + SA{x) 

B{x) B{x) + 5B{x) . (14) 

The dynamics of the disturbances 5 A and 5B is deter- 
mined by the equations of motion (7). In a general local 
inertial frame, we obtain the following linearized equa- 
tions: 

= U5A - 2B^'df,5A + 2m\Bl8A 

- 2At'd^6B + (2A^S^ + 4m\AbBb)6B , 

= DSB - 2Bi'd^5B + 2Ai'd^5A 

+ (6m|s2 - A^'A^ + B''B^)SB . (15) 



Here, the cosmological axion-dilaton background is de- 
scribed by Ab = A{x), Bb = B{x), A'' = Bj^^S''^, 
and B^' = B^^df'B. Note that Ay, and B^ satisfy the 
equations of motion (7). Locally, the above definitions 
yield effectively constant quantities because the axion- 
dilaton background varies appreciably only on cosmolog- 
ical scales. In what follows, we can therefore take A^, 
Bb, Af^, and B^ as evaluated at x = a;o. 

An ansatz with plane waves exp(— ip • x) gives a set of 
algebraic momentum-space equations, as ususal. These 
two equations are governed by the characteristic matrix 
C: 

^_ f p^-M%-2iB.p -2iA-p-2M^ \ . . 

2iA-p P^-Ml- 2iB-p J ' ^^^> 

where Ml = 2m\Bl, Ml = 6m%Bl - A- A + B B, and 
= A-B + 2m\AbBb. With the input values (10), 
and are positive in the recent past of our model 
universe. However, other parameter choices can lead to 
negative values for M^ and M^. A detailed discussion of 
such a situation would be interesting but lies beyond the 
scope of the present work. 

The plane-wave dispersion relation describing the 
propagation of the various modes is determined by 
det(C) = 0, as usual. We obtain 

= {p^ -Ml- 2iB-p){p'^ -Ml- 2iB-p) 

+ A{iM'^ -A-p) A-p . (17) 

The imaginary terms in this dispersion relation are a 
direct consequence of the non-Hermiticity of C. They 
will lead to exponentially increasing or decaying solu- 
tions. This feature is characteristic in cases with broken 
spacetime-translation symmetry and the associated non- 
conservation of 4-momentum. 

Another important feature of the dispersion relation 
(17) are terms in which the plane- wave 4-momentum p^ 
is contracted with A'^ and B'*. Since these 4- vectors 
are taken as a constant nondynamical background in 
the present context, they select a direction in local in- 
ertial frames violating particle Lorentz invariance. As 
in the case of varying couplings, this result is intuitively 
reasonable: the effective vacuum containing the axion- 
dilaton background acts as a spacetime- varying medium 
in which the disturbances 6 A and dB propagate. More- 
over, the spacetimc-dependent cosmological solution for 
the scalars breaks translation symmetry, while transla- 
tions and Lorentz transformations are intertwined in the 
Poincare group. The above argument for Lorentz break- 
ing is not only confined to the present supergravity cos- 
mology. It applies to any model in which scalar fields ac- 
quire expectation values varying on cosmological scales: 
the propagation of localized disturbances in the back- 
ground expectation value will typically violate Lorentz 
symmetry. Note that scalars with a "rolling" time de- 
pendence, such as those in infiation, quintessence, and 
k-essence models, play a key role in cosmology. 
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The tight experimental contraints on the modification 
parameters in the SME [35] imply that possible Lorentz 
and CPT violations must be minuscule for known par- 
ticles. For instance, the best bounds for radiation and 
matter can be found in Refs. [36,37] and [38], respec- 
tively. In the case of Lorentz-breaking dispersion rela- 
tions, kinematical threshold analyses are a useful tool 
[39]. Investigations of ultrahigh-energy cosmic rays, for 
example, can place Planck-scalc bounds on some SME 
coefficients for Lorentz violation [40]. The above exper- 
imental methods are in principle also applicable in the 
present context. However, in the absence of direct ob- 
servational evidence for such fundamental scalars, the 
present type of Lorentz violation is phenomenologically 
less interesting at the present time. Moreover, the dis- 
cussion in the previous section implies an expected size 

0{Ho) for the background gradients, which is below 
current experimental sensitivities. 

Next, we take advantage of coordinate Lorentz invari- 
ance and continue our analysis in a local comoving in- 
ertial frame. This is also appropriate from an experi- 
mental viewpoint because in realistic situations laborato- 
ries move nonrelativistically with respect to such frames. 
Then, the scalars A and B depend only the time t [41], 
so that = {A, 0) and = {B,Q), and the dispersion 
relation (17) can be solved for = p^p^: 

p^ = ^{Ml + M%) + 2iBE 

± \l\{Ml - MlY - 4AE{iM-^ - Ae) . (18) 

Here, we have denoted p'^ = E. 

In principle, the dispersion relation (18) can be cast 
into the conventional form E{p), which involves in the 
present case the general roots of a quartic equation. How- 
ever, the expressions for the plane- wave energies become 
more transparent in certain limits. For example, we can 
consider a case in which the non-standard dispersion- 
relation terms are small, i.e., 

A^A^', B^B'', Af2 « Ml, M| . (19) 

Although the two scales are intertwined in the equa- 
tions of motion (7), this situation can in principle be 
realized. Then, in the nonrelativistic limit ]p] ~ 0, 
depends only on the above two scales, and one expects 
\E\ < 0{Ma, Mb, ). It follows that the dispersion rela- 
tion (18) reduces in leading order to 

p' = i(Ml + M|) ± i(Ml - Ml) . (20) 

In such a situation. Ma and Mb would therefore be the 
effective axion and dilaton masses. We remark that the 
condition (19) is violated in the model with the param- 
eters (10). This case is discussed below when positivity 
and causality are investigated. 

Another useful approximation is the ultrarelativistic 
limit 



E+ ~ ±\p\+A„+iBo , 

E_ c:^ ±\p\ - Ao + iBo . (21) 



Here, the subscripts -|- and — correspond to signs of 
the square root in Eq. (18). Asymptotically, the real 
part of the energy variable grows linearly with the 3- 
momentum, as in the conventional case. Note, how- 
ever, the presence of imaginary terms. Depending on the 
sign of Bq, they lead to exponentially growing or decay- 
ing solutions implying non-conservation of 4-momentum. 
This feature does not come as a surprise because the 
cosmological axion-dilaton background breaks spacetime- 
translation invariance. One might argue against expo- 
nentially growing solutions as being unphysical. How- 
ever, it is important to observe that lm{E) is determined 
by the minuscule Lorentz-violating terms and does not 
increase with the 3-momentum. Moreover, our solutions 
are valid only within a localized spacetime region with 
a characteristic size At satisfying BAt ^ 1. Since the 
local time r obeys At < At, the change in amplitude de- 
scribed by exp(SAT) ^ 1, and thus the resulting effect, 
must be small. One might also suspect that the imagi- 
nary contributions to the energy can lead to interpreta- 
tional difficulty when treated within quantum field theory 
(QFT). For instance, the free QFT Hamiltonian appears 
to be non-Hermitian and Feynman boundary conditions 
for the Green function may seemingly no longer be im- 
posed. However, these apparent obstacles can be avoided 
if the time-dependent background is treated perturba- 
tively. Then, particle number, for instance, fails to be 
conserved, which is consistent with the above plane- wave 
solution of varying amplitude. 

Lorentz-breaking dispersion relations typically lead to 
positivity problems, super luminal group velocities, or 
both at high energies [32,42]. In the present model, the 
situation is further complicated by the presence of imagi- 
nary terms, so that we restrict ourselves to a few remarks. 
As before, we work in a local comoving inertial frame, so 
that rotational symmetry is maintained. 

Figure 7 depicts the real part of the energy variable 

clS Si function of a momentum component in a model 
with parameters (10). The dotted line represents the 
momentum-space lightcone. As expected, there are four 
branches, each corresponding to one of the four roots of 
the dispersion relation (17). Note that in the plotted mo- 
mentum range two of these branches, shown as dashed 
lines, lie in the shaded region of spacelike momenta. The 
presence of spacelike 4-momenta is normally associated 
with negative particle energies in some inertial frames. 
However, instabilities are absent in the present model: 
the conserved axion-dilaton energy including the back- 
ground remains always real and positive definite, as is 
evident from the energy-momentum tensor appearing in 
Eq. (6). 



7 




-2 -1 ft 1 2 



FIG. 7. Real part of the energy E versus a momentum 
component for the model parameters (10). The dotted lines 
correspond to the momentum-space lightcone. The solid and 
dashed branches are associated with the various roots of the 
dispersion relation (17). Note the presence of spacelike mo- 
menta and superluminal group velocities. 
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FIG. 8. Imaginary part of the energy E versus a momen- 
tum component for the model parameters (10). The dotted 
lines represent the momentum-space lightcone. The solid and 
dashed lines correspond to the respective solid and dashed 
branches in Fig. 7. For large values of the momentum, the 
solid and dashed lines lie on top of each other. The region of 
spacelike momenta is shaded. 

Asymptotically, the magnitude of the group velocity 

Vg = Vp- E{p) equals unity because the conventional 
Lorentz-symmetric term {p^p^i)'^ dominates the disper- 
sion relation (17) at high energies. However, inspection 
of the slope of a dashed branch in Fig. 7 reveals the 
possibility that |iTg| > c at low energies. For a back- 
ground that is fixed and fundamentally nondynamical, 
such superluminal group velocities violate microcausal- 
ity [32]. But for conventional dynamical backgrounds 
leading to group velocities exceeding c, causality is typi- 
cally maintained. Examples in ordinary electrodynamics 
can be readily identified: the occurrence of superluminal 
Ug in amplifying media with inverted atomic populations 
has been verified theoretically and experimentally [43], 
but causality violations are absent [44]. The underlying 



Lorentz-covariant dynamics of the our supergravity cos- 
mology leads therefore to the conjecture that causality is 
maintained in the present context. The apparent super- 
luminal particle propagation would then be an artifact of 
the approximations employed in the derivation of the dis- 
persion relation (17). A complete investigation of these 
issues would be interesting but lies well outside the scope 
of this work. 

For completeness, the imaginary part of the energy 
variable as a function of momentum is plotted in Fig. 
8. The momentum-space lightcone is again represented 
by the dotted line. The shaded region is associated with 
spacelike 4-momenta. The solid and dashed lines corre- 
spond to the solid and dashed branches in Fig. (7). At 
large 3-momenta, all four branches overlap and Im(S) is 
bounded, which is in agreement with Eq. (21). 



VI. SUMMARY 

This work has considered cosmologies with scalar fields 
that are motivated in candidate fundamental theories. In 
such a context, the scalars typically acquire varying ex- 
pectation values that can be associated with an acceler- 
ated expansion of the universe, varying couplings, and 
Lorentz- violation. 

More specifically, we have investigated an TV = 4 su- 
pergravity model in four dimensions. In this framework, 
standard plausible arguments lead to a variation of the 
axion and the dilaton on cosmological scales. As a result, 
the propagation of the axion and the dilaton is governed 
by a Lorentz-violating effective dispersion relation. We 
expect this feature to be generic in models with scalars 
varying on cosmological scales. 

The axion-dilaton background also affects the expan- 
sion of the universe and results in spacetime-dependent 
electromagnetic couplings. Model parameters exist that 
lead to a behavior of the scale factor consistent with the 
observed late-time cosmological expansion. The variation 
of a implied by this parameter set lies mostly outside ex- 
perimental constraints. However, the time dependence of 
a is roughly in the order of magnitude suggested by the 
Webb data and displays desirable nonlinear features. 
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APPENDIX A: ADDITIONAL REMARKS 
REGARDING THE SPECTRUM OF SOLUTIONS 

Wc begin our considerations with the observation that 
for all nontrivial physical input values the universe is 
always expanding in our supergravity cosmology. Sup- 
pose this were not the case. Then, at some time t = to, 
there must be a transition from the presently observed 
cosmological expansion to a subsequent epoch involving 
contraction. This implies d{to) = at the time to by con- 
tinuity of d{t). Then, the first of the Einstein equations 
(6) yields at t = to 



2B2 



+ 2p = Q 



(Al) 



which can only be satisfied in the trivial situation of an 
empty universe with zero total energy density or in un- 
physical cases involving negative matter densities p < 0. 

We continue our discussion by considering special 
choices for the parameters uia and tob- The situation 
rriA = and ms = has been discussed previously in 
the literature [16]. Let us therefore focus first on the 
case niA ^ and tob = 0. With this choice, we find the 
following analytical solution: 



a{t) = M^t^'^ , 
A{t) = ±{mAt)-^ , 
B{t) = ±(mAi)"' • 



(A2) 



Here, the signs in the equations for A{t) and B{t) can 
be selected independently. The absence of integration 
constants indicates that more general solutions exist (cf. 
model P2 below). Note that the cosmological expan- 
sion implied by Eqs. (A2) is essentially that of a conven- 
tional matter-dominated universe, which always deceler- 
ates. The time dependence of the fine-structure parame- 
ter is given by 



mit^ 



1 



An 2mAt + ■m\t^ 



(A3) 



It follows that at late times, a increases linearly with t. 

Next, we look at the case niA = and ms 7^ 0. One 
can verify that 



a{t)= ylcnt^^ 
A{t) = Ao + Ast-^ , 
B{t) = ±\/2(msi)"^ 



(A4) 



is an asymptotic {t — > 00) solution. The integration con- 
stants Aq and A3 can be chosen freely. This sohition also 
describes a non-accelerated cosmological expansion. For 
Ao ^ ±1, the late-time behavior of a{t) is given by 



a{t) = 



1 {1-Alf 



4:V2tt 1 + Al 



mst 



(A5) 



In a situation with Ao = ±1, the fine-structure parame- 
ter decreases asymptotically as follows: 



(A6) 



We also mention that with the choice A3 = 0, Eqs. (A4) 
become an exact solution. In this case, the time-evolution 
of the electromagnetic coupling reads 



ait) = 



1_ 4 + 4(1 + ADmlt'' + (1 - Al)^m%t^ 
47r l^ruBt \/2(l -h Ag)m|i3 



(A7) 



The general situation, in which both and ms are 

nonzero, seems to evade systematic analysis with analyt- 
ical methods. Numerical investigations reveal the exis- 
tence of a broad spectrum of qualitatively different so- 
lutions. Figure 9 conveys a flavor of the diversity in ac- 
celeration behavior of our supergravity cosmology. The 
second derivative of the scale factor, a(t), is plotted ver- 
sus the comoving time t for the four parameter sets PI to 
P4 given in Table I. The shaded region marks accelerated 
cosmological expansion. The input values PI arc associ- 
ated with a cosmology of permanent deceleration, qual- 
itatively analogous to a conventional matter-dominated 
universe without cosmological constant. The model pa- 
rameters P2 lead to an initial deceleration period followed 
by an asymptotic accelerated expansion. This behavior 
agrees qualitatively with that of the phenomenological 
model (10). An overall decelerated expansion with a 
transient period of acceleration results from the input 
P3. The parameter set P4 implies an initial deceleration, 
followed by transient periods of acceleration and deceler- 
ation before the expansion continues to accelerate. 

It is worth emphasizing that some of the solutions we 
have obtained exhibit the interesting feature of a tran- 
sient accelerated expansion (cf. model P3). This is a pos- 
sible fix for the recently discussed problem concerning an 
eternally accelerated cosmological expansion. Such uni- 
verses would pose a challenge for string theory, at least in 
its present formulation, since string asymptotic states are 
inconsistent with spacetimes that possess event horizons 
in the future [45]. 




" b 



12 3 4 

comoving time ( in 10"*^ GeV' 

FIG. 9. Acceleration a{t) of the scale factor versus comov- 
ing time t for the various input values given in Tabic I. The 
shaded region corresponds to accelerated expansion. The di- 
versity in qualitative behavior is apparent. 
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TABLE L luput-parameter sots PI. P2. P3. and PI. 



Parameter 


PI 


P2 


P3 


P4 


ruA in 10-^2 Q^y 





1.5 





1 


TJiB in lO-'*^ GeV 


10 





100 


100 


c„ in 10-84 GeV^ 


2 


2 


2 


2 


a (in) 


1 


1 


1 


1 




1.023 


1.023 


1.023 


1.023 


i(i„) in 10-47 GcV 


47 


47 


47 


-100 


B(in) 


0.022 


0.022 


0.022 


0.022 


B(tn) in 10-45 GeV 


-25 


-25 


-25 


-60 


tn in 1040 GeV-^ 


56 


51 


54 


51 
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